THE FAA DI BRUNO FORMULA REVISITED 
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Abstract. We present an intuitive approach to (a variant of) the Faa di Bruno formula which 
shows how this formula may have been (re)discovered many times in history. Our representation 
' for the n-th derivative of the composition / o g of two smooth functions / and g on R uses a simpler 

summation order so that the mysterious condition 61 + 2&2 + ■ ■ ■ + nb„ = n in Faa di Bruno's formula 
does not appear. 
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1. HOW ONE MAY (RE)DISCOVER ONESELF THIS FORMULA 

Let / and g be two functions on K. whose n-th derivatives exist. The first 4 derivatives of / o g 
are easy to calculate: 



(fog)' = (f'og)g' 

(/ojpl = (/ (2) og)( 3 ') 2 + (/'o 9 )g< 2 ' 

(/o 9 )< 31 = (/ (3) o S )( S ') 3 + 3(/< 2 >o 9)9 < 2 y + (f'og)g^ 

(/°9) (4) = (/ <4) og)(g') 4 + 6(/< 3 > o g )g( 2 >( 9 ') 2 + 4(/' 2 > o g)c 

£f~) • The expressions get rapidly longer and more complicated. For example we have 42 summands 

for n = 10. Note that the number of terms equals the partition number p(n), that is the number 
of ways (without order) of writing the integer n as a sum of strictly positive integers; by the 

Hardy-Ramanujan formula we have p(n) ~ ^J^- e 77 "^ 2 "/ 3 . 
Let 

M j = {k = (k u . . . , kj) e h>k 2 >--->kj>l} 
be the set of ordered multi-indeces in N* = {1, 2, . . . }. If g is a function defined on R, and if is 
the n-th derivative of g, then we denote by g^ the function Y\{=i 9^ ki \ where k = (k±, . . . ,kj) € 
M.K Also, is, by convention, equal to the function g itself. 
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The classical Faa di Bruno formula from ca. 1850 gives an explicit formula for (/ o g)^: 

where the sum is taken over all different solutions in nonnegative integers b\ , hi , . . . , b n of 

bi + 2b 2 H h nb n = n and j := b\ H V b n . 

A nice historical survey on this appeared in [2]. See also [T]. 

Without being aware of that formula, I developed around 1976-1980 the following formula: 

(1-1) (/°5) (n) (*) = E/ ( %(*))( £ C%gM(x)\ 

i =1 k e M j 

\k\=n 
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where 

' n 



k 



~[N(k,', 



■3 ~ n > 



(ft \ Tl ^ 
, I = — 71—, r~r) where \k\ := k\+- • -+k 
k) fci'A^' • • • kj\ 

and N(k, i) is the number of times the integer i appears in the j-tuple k (i £ N* and k £ (N*) J ). 

For example, the coefficient C^ 4 1 ^ of the term g^(g') 2 when looking at the 6-th order derivative 
of / o g is 

C 6 - 1 6! - 15 
^(4,1,1) - 2 ! 4! . l! • 1! ~ 

and the coefficient C^ 2 1 1 1 ^ of the term g^g 1 ' '(<?') 4 i n the 10-th derivative is 

C 10 - I 101 _ 3150 

u (4,2,i,i,i,i) ~ 4! 4! . 2! • 1! • 1! • 1! • 1! ~ 

The difference between our formula and the Faa di Buno formula is that we use a simpler 
summation order and do not consider exponents of the form &-• = 0. In particular, we do not need 
summation over those (&i, . . . , b n ) satisfying (the difficult to grasp) condition Yl7=i = n - That 
these two formulas are equivalent though, immediately follows from a direct comparison of the 
coefficients. Indeed, for fixed j and k = (k±, . . . ,kj) E |fc| = n, we have: 

n\ \ 1 1 / , 1 1 \ 1 

n! 



~ \h\---kAj N(k,l)\ N(k,n)l ~ »/!•••*<!,/ 



6j. times 



n\ rr / 1 



n 



where the bi m are those exponents that are different from zero and where k has been represented 
in the canonical form k = . . . , i\, ■ ■ ■ . . . , ij ) in decreasing order. Note that 

bi 1 times 6^ times 

n t j 

i=l s=l p=l 

and that 



^ &i = bn H 1- &i £ = i- 



Next I would like to present the (intuitive) steps that led me to the discovery of the formula 
above, at pre-PC times; the first (non-programmable) slide rule calculator SR50 had just appeared. 

1) I calculated explicitely the derivatives (/ o g)^ up to the order 10 and wrote them down in 
a careful chosen order (see figure 1); 

2) An immediate guess is that 
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(/o0)W(x) = £/Wfc(x)) E 4^ (k) ( 



i =1 k e M j 

\k\=n 

for some coefficients to be determined. 

3) Next I gave an inductive proof that this representation is correct; that needs the main step 
of the construction: where does the factor • • • with fa + • • • + kj = n + 1 comes from? So 
let us look at the ordered j-tuple (fa,... , kj). This tuple is generated, through anti-differentiation, 
by the j j-tuples 



that is, 



(fa - l,k 2 , ■ ■ -,kj), (ki,k 2 - l,fc 3 ,... , kj), , (fa, ...,kj- 1); 

g{kl-l) g {k 2 ) . . . g {k 3 )^' = g(kl)g(k 2 )... g (k j ) + ... j 

5 (fci) 5 (fca-l)... 5 (fci)y = g (ki) g (k2) ... g (kj) + 
etc. 

The main difficulty being that the components kj are not pairwise distinct. So their "multiplic- 
ities" had to be taken into account. This lead to the guess that one may have 

c\ = l 

+1 = A N(k-4,ki-l) 

k N{k ,h) k-4 

where for i = ej = (0, . . . , 0,^,, 0, . . . , 0), e{ G N j , k e M j ,\k\ = n+1 and 1 < j < n+1. 

i—th 

Note that the j-tuple k — eP- is not necessarily represented in the canonical form with decreasing 
coordinates. Also, if the i-th coordinate of k is one, then the z-th coordinate of k — is and we 
identify k — e J { with the associated (j — l)-tuple. 
For example in the case k = (3, 1, 1, 1) we have 

6 1 5 1 5 1 5 ^ 5 

20 = C( 3]ljll ) = j-C(2,i,i,i) + 3^(3,0,1,1) + 3^(3,1,0,1) + 3^(3,1,1,0)' 

where (3, 0, 1, 1), (3, 1, 0, 1) and (3, 1, 1, 0) are identified with (3, 1, 1). 
With these recursion formula the inductive proof went through. 

4) Next one has to guess the explicite value of c^, \k\ = n. Now there are (£) = ways to 
choose k\ objects out of n, then k 2 objects of the remaining ones, and so on. Due to the multiplicty, 
one has again to divide by N(k,i)\. 

CD 

This gives the guess that = v 7 . 

k l[N(k,i)\ 

i 

5) These coefficients actually satisfy the recursion relation above. Since c\ = 1, and the fact 
that the recursion relation determines uniquely the next coefficients, we are done: = c^. 
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2. Further formulas and questions 
Applying formula for the function /(x) = logx,x > and g(x) = e x gives 



C£=0, n>2 ; 

i =1 k e M j 

\k\=n 

whereas for /(x) = x n and <?(x) = e x one obtains 



e 

|fc|=n 

In particular, L := C£ < n n . Is there an explicit expression for LI If one uses f(x) = g{x) = 

k:\k\=n 

e x , then 

(e^) {n) \ x=0 = eL. 
One may also ask the following questions: 

(1) What is C k 0-<3< n V 

k £ M j 
\k\=n 

(2) What is max{C£ : \k\ = n}? 

(3) Is there a formula for the number of partitions of n with fixed length j? 

In our scheme (figure 1), one can give easy formulas for the coefficients in each column. In fact, 
each element in a fixed column is a multiple of the first coefficient. More precisely, if k\ > k-i > 
kj > 1, then we have: 

^tl.^l, ••-,!) =<7 (V.,*y)( \ )' 

t — times 

We observe that several columns coincide; for example C^ 2 2 2 ) = ^{i 2) an< ^ so ^ e elements of 



the associated columns are the same. 

vc nf r*r\ I n t~y~> c f nof nmnriirlo I nncf nco + ri o +- # 

(2,2,2,i) (4,2,i) 



There are actually infinitely many pairs of colums that coincide (just use that C^oo 1 = 



for every i > 5. 
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Are there triples (or higher number) of columns that coincide? 

Acknowledgements I thank Jerome Noel for his help in finding the exact T^X-commands for 
creating the tabular form above. I also thank Claude Jung for valuable discussions around 1980. 
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